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ABSTRACT 

We present the angular power spectrum of the CMB component extracted with Fas- 
tICA from the Background Emission Anisotropy Scanning Telescope (BEAST) data. 
BEAST is a 2.2 meter off-axis telescope with a focal plane comprising 8 elements at 
Q (38-45 GHz) and Ka (26-36 GHz) bands. It operates from the UC White Mountain 
Research Station at an altitude of 3800 meters. The BEAST CMB angular power 
spectrum has been already calculated by O'Dwyer et al. using only the Q band data. 
With two input channels FastICA returns two possible independent components. We 
found that one of these two has an unphysical spectral behaviour while the other is a 
reasonable CMB component. After a detailed calibration procedure based on Monte- 
Carlo (MC) simulations we extracted the angular power spectrum for the identified 
CMB component and found a very good agreement with the already published BEAST 
CMB angular power spectrum and with the WMAP data. 

Key words: methods - data analysis - techniques: image processing - cosmic mi- 
crowave background. 



1 INTRODUCTION 

The recent outstandinK results from WMAP satellite 
iBennett et a l. 2003) have definitely put us into the era of 
precision cosmology with an accurate determination of the 
CMB angular powe r spectrum up to I ~ 800. In addition the 
DASI experiment jLeitch et al.ll200^ has clearly reported 
detection of E mode CMB polarization. The situation will 
improve even further with the Planck satellite, a third gen- 
eration CMB space mission which will map microwave emis- 
sion over the whole sky with an unprecedented combination 
of angular resolution and sensitivity. In the meantime a ple- 
tora of both ground-based and balloon borne experiments 
will produce accurate measurements (better than WMAP) 
over limited sky regions. 

Today one of the main limitations to the accuracy is 
the presence of other astrophysical sources between us and 
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the Last Scattering Surface, which contribute to the mea- 
sured signal. These foreground contaminants consist mainly 
of Galactic emission (synchrotron, free-free and dust emis- 
sion), compact galactic and extragalactic sources, and the 
Sunyaev-Zel'dovich effect from cluster of galaxies. The chal- 
lenge is to identify and remove such foreground emissions 
with high accuracy and reliability in order to obtain cleaned 
CMB maps. This is crucial in deriving precise cosmological 
information from the CMB power spectrum. 

Many works have been dedicated to component 
separation, and different algorithms have been pro- 
posed. Traditional separation t echni q ues, from Wiener fil- 
tering jTe gmark fc Efstathiouj Il996l iBouchet et"al] 119991 : 
iPrunet et al. 2 001) t o Maximum Entropy Method (MEM) 
jHobson et allligga: IStolvarov et~ai]|2002h . have generally 
been employed. They achieve good results, but they re- 
quire prior knowledge about the signals to be separated (e.g. 
spatial templates and frequency dependance), whereas the 
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avaible full-sky foreground priors actually are not completely 
reliable. 

Recently a blind separation approach has been devel- 
oped, which works without the need of priors, except about 
the statistical features of the components. Indeed this tech- 
nique, based o n the Independent Component Analysis (ICA) 
^Comonlll994^ . exploits the statistical independence of the 
sky siRuals. It wa s first implemented as a neural network 
iBaccigalupi et al.n2000^'l , and then optimized in a fast algo- 
rithm. FASTlCA~ lMaino et al which was successfully 

tested on simulated sky maps similar to those that Planck 
will produce. FastICA has shown good perf ormance also 
when applied to real data from CO-B-E-DMR jMaino et alJ 
|20o3), with results on CMB anisotropy and foreground con- 
tamination consistent with previous and independent anal- 
yses. 

In this work we apply FastIC A to another real data set, 
from the Background Emission Anisotropy Scanning Tele- 
scope (BEAST). In SectionHwe briefly recaU the main fea- 
tures of the FastICA approach and its assumptions. In Sec- 
tion|31 after describing the BEAST instrument and the maps 
produced, we explain the procedure followed to apply Fas- 
tICA to BEAST data. The results obtained are presented 
in Section [3.11 and the CMB reconstruction quality, tested 
with Monte Carlo simulations, is analysed in Section 13.21 
Sections |1| deals with the normalisation of the CMB sig- 
nal extracted with FastICA. In Section |S] we extract the 
FastICA CMB spectrum. Finally, a critical discussion and 
conclusions are presented in Section |S| 



2 COMPONENT SEPARATION WITH FastICA 

Before describing the FastICA method, it is useful to recall 
briefly the data model FastICA refers to and the principal 
assumptions from which it derives, described in detail in 
iMaino et al"! (|20o3). 

Let us suppose that the sky radiation, as a function of 
direction r and frequency u, is a superposition of A'^ different 
signals Sj[r,v) and that it is observed by an experiment 
with M frequency channels whose beam pattern is B[r, v). 
Let us further suppose that, for each signal, frequency and 
spatial dependence can be factored into two separated terms, 
fiiv) and Sj[r) respectively, and that B is shift-invariant 
and frequency-independent. Then the data model can be 
written as: 

x(r) = As(r) * B{r) + e(r) = As(r) -I- e(r) , (1) 

where each component, Sj, of the vector s is the correspond- 
ing source function convolved with the B beam pattern. The 
matrix A is the mixing matrix, which includes the frequency 
response, and e(r) is the instrumental noise term. 

The FastICA algorithm obtains both the mixing ma- 
trix A and the signals s from observed data x assuming 
that 

• the signals s are independent random processes on the 
map domain; 

• all the signal, but at most one, have non-Gaussian dis- 
tribution. 



application i n an astrophysical context is described in 
iMaino et al.l (l2002l'l . Independent components are extracted 
maximising a suitable measure of non-Gaussianity that is 
robust against noise: this is the the so-called neg-entropy. 

FastICA estimates separation matrix W row by row, 
maximizing the non-Gaussianity of the component w-^x, 
where is a row of W, such that the trasformed vari- 
ables y = Wx are the independent components. 

In particular the FastICA algorithm operate s with 
a neg-entropy approximation armen fcOkl 1200(1 : 

iHvvarinenI 11999) . which can assume three different forms, 
depending on the regular non-quadratic function chosen in 
its expression: g(u) = v?' , g{u) = u exp(— u^) or g{u) = 
tanh(?i), where u — w^x. In the following we indicate these 
functions as p, g and t respectively. The best choice of the 
function depends on the statistical proprieties of the com- 
ponents: kurtosis, or p, can be used for sub-Gaussian com- 
ponents in absence of outliers; g may be better when the 
components are higly super-Gaussian or when robustness i s 
important; t is a general-purpose function llHvv arme nll999ll . 
However we do not know a priori the statistics of the inde- 
pendent signals. 

Once the separation matrix W is obtained, since we 
have X — W~''"y, we can derive the frequency scalings for 
each independent component: the scaling between u and v' 
of the j**^ component is given by the ratio of /W~,\ If 
the spectral behaviour is given by a power law with index 
P, then we have 13 = \og[W~^ /W~]]/ \og{v /v'). Therefore a 
negative value of frequency scaling indicates that the recov- 
ered component does not have a physical behaviour. 

Furthermore we can estimate the noise in the recon- 
structed maps. If we perform noise constrained realizations 
Hx for each frequency channel, the corresponding noise re- 
alizations in the FastICA outputs are given by Wnx. 



3 APPLICATION TO BEAST DATA 

The Background Emission Anisotropy Scanning Telescope 
BEAST) is a 2.2 m eter off-axis Gregorian telescope 
Figueiredo et al .120051) with a focal plane consisting of 6 Q- 
band (38-45 GHz) and 2 Ka-band (26-36 GHz) corrugated 
scalar feed horns cou pled with cryogenic HEMT ampliflers 
JChilders et alJl2005ll . The instrument was installed at the 
UC White Mountain Research Station at an altitude of 3.8 
km in 2001 July. Data considered here come from two dif- 
ferent campaigns: one until December 2001 and the second 
one in 2002 (February and August /September). 

BEAST produced two maps covering an annular sky re- 
gion around the NCP from 33° < 5 < 42° with a resolution 
of 23' in Q-band and 30' in Ka-band. The s ky maps are 
pixelized according to the HEALPix^ scheme l|G6rski et alJ 
1999) with a resolution parameter Ngide — 512 correspond- 
ing to pixel size of 6.9'. 

For a proper analysis of the results that FastICA will 
obtain, we need to know the instrumental noise properties, 
namely noise level in the two frequency bands a nd its spatial 
distribution. To estimate instrumental noise uMeinhold et alJ 
I2OO5I have made "difference" maps at the two frequency 



A detailed explanation of th is strategy can b e found in 
iHvvarinen fc Oial Jl997ft and iHvvarinenI (Il999h . while its 
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Table 1. ICA frequency scaling of the astrophysical component from BEAST maps smoothed at the three angular resolution, and for 
the three galactic cuts. 
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Figure 1. left: The BEAST Q-band map smoothed to 40'; right: 
galactic plane is removed for \b\ ^ 17.5°. 



bands separately. They binned data from first half of obser- 
vation into one map and the second half into another one. 
For each band the "difference" map is the pixel by pixel dif- 
ference of these two maps, divided by 2 to maintain noise 
statistics as in the sum map. Therefore these maps should 
not contain in principle, any sky signal, but only noise and 
the root mean square (rms) of these maps is a measure of 
instrumental noise. The S/N ratio is quite poor: for the Q- 
band map Moinhold ct al. 2005 have found a value ~ 0.11 
at 23' resolution, which becomes ~ 0.57 when the map is 
smoothed at 30'. The Ka-band map shows a higher noise 
contribution. In addition from the "difference" maps we can 
see that noise is gaussian but, due to the scanning strategy, 
is not uniformly distributed on the sky. 

Before applying FastICA a smoothing of the maps 
is required, due to the different angular resolution in the 
two bands, since in FastICA approach it is assumed a 
frequency-independent beam pattern (see Eq. . We have 
smoothed the Ka and the Q maps to the same angular reso- 
lution, choosing values of 30', 40' and 60' in order to increase 
the S/N ratio. To obtain maps with significant signal level 
we decided not to smooth to resolution greater than 60', be- 
cause of the 10 Hz high-pass filter applied to the BEAST 
data in the reduction processing: this indeed produces a sig- 
nal cutoff on angular scales ^ 6°. 

We applied FastICA to the BEAST maps at the three 
resolutions working with all the non-quadratic functions de- 
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The reconstructed astrophysical component (with g function). The 



scribed before identified by p, g and t respectively. With 
two input sky maps FastICA is able to reconstruct only 
two outputs, since the estimated matrix W has dimension 
2x2. For all the considered cases one of the two is an as- 
trophysical component, while the second component, due to 
the low S/N ratio, always has a negative frequency scal- 
ing, which suggests unphysical behaviour as explained in 
the previous section. Performing a component separation 
on the full observed sky FastICA recovers an astrophysical 
signal with a frequency scaling consistent with free- free emis- 
sion. This is due to the strong galactic emission in the two 
plane crossings. Therefore, in order to reconstruct the CMB 
component, we cut the sky regions where foreground emis- 
sion overcomes CMB. We remove the three strongest point 
sources and the galactic plane, cutting out data within three 
different valu es of galactic lat itude: |&| ^ 17.5°, |&| ^ 20° and 
|6| ^ 22°. Mei ia et al.l ^2005^ have estimated that, removing 
from BEAST maps regions with |6| 5C 17.5°, the individual 
galactic contributions remain below ~ 1% of the map rms. 



3.1 The CMB component 

We wanted to verify if the component extracted by FastICA 
is indeed consistent with CMB signal and the first figure of 
merit is the expected frequency scaling of the CMB between 
the two BEAST frequencies. BEAST data are in antenna 



4 S. Donzelli et al. 



30 arcminutes 



40 arcminutes 



0.70 


- 




0.70 






0.65 






0,65 


- 


























0,60 






0,60 


















0.55 






0.55 


















0.50 




0.50 




* 











60 arcminutes 



','\ 

h : 4 




- 























0.0 0.5 1.0 1.5 2,0 2.5 3.0 5.5 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 
frequency scoling 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 



Figure 2. Correlation between ICA CMB and simulated CMB in Ka (plus sign = p, square = g, 'x' = t) and Q band (asterisk = p, 
diamond = g, triangle = t) towards recovered frequency scaling at 30', 40' and 60'. Line-dot lines show the expected scaling (1.022). Dot 
lines indicate the minimun correlation and the maximum scaling for a good CMB reconstruction. We can also observe that correlation 
with CMB-Q map is greater than with CMB-Ka map. 



temperature and CMB fluctuations STa.cms are related to 
brightness temperature fluctuation ^Tcmb by: 

5Ta,cm^{v) = _ (JTcMB , (2) 

where x = hu/kTcMB- Assuming for the CMB a bl ackbody 
temperature Tcmb = 2.725 K ( Mather et al,lll999l) . we ex- 
pect for the CMB component a frequency scaling between 
30 (Ka) and 41,5 (Q) GHz equal to 1,022, 

After removing the sky regions strongly contaminated 
by galactic emission, FastICA recovers an astrophysical 
component with frequency scaling quite different from that 
expected for CMB, as shown in Table Q This is due to the 
high noise contribution, which aflects also the recontructed 
astrophysical signal. Furthermore from Tabled we can ob- 
serve that results obtained with t function are in general 
worse than p and g results. Finally, increasing angular scale, 
with p and g the frequency scaling approaches the expected 
value in most of the cases. Nevertheless the spatial pattern 
of the recontructed astrophysical component resembles the 
Q map pattern, as shown in Fig, for the 40' smoothing 
case. 

After this first indication we proceeded by verifying 
that reconstructed maps at different resolution are consis- 
tent one each other. After smoothing all maps down to 60', 
we have calculated the correlation between maps with dif- 
ferent original resolution, finding Spearman correlation co- 
efficients rs > 0,8 (with the exception of correlation between 
maps 30' — 60' obtained with t function, for which Vs ~ 0,7), 
This correlation indicates that FastICA recovers the same 
astrophysical signal at every resolution. Furthermore we ob- 
served that all the reconstructed maps have high spatial 
correlation with the Q band map smoothed at the same res- 
olution, with Ts > 0,9 except for t results at 40' {r, ~ 0,7) 
and at 60' (r^ ~ 0,6), This is indeed expected since the S/N 
ratio is larger in Q band than in Ka, Finally, we verified 
that there are not significant changes in the astrophysical 



component reconstruction when extending the galactic cut. 
Indeed reconstructed maps applying different cuts are con- 
sistent one each other (r^ ~ 0,9), Also the frequency scal- 
ing does not change significantly with the galactic cut (see 
Table and there is not a well defined trend in the scal- 
ing variations with the cut extension. This points out that 
the galactic contribution in the recontructed astrophysical 
component is not relevant. Therefore we are confident that 
FastICA recovers a signal dominated by CMB anisotropics. 



3.2 Testing results with Monte Carlo simulations 

Given the poor S/N ratio in the BEAST data, we prefer to 
test the performance of FastICA by simulation. In order 
to analyse CMB reconstruction quality we performed 100 
Monte Carlo simulations in which sky signal is simulated, 
observed following BEAST observing strategy, reduced as 
the actual data and then analysed by FastICA, We already 
demonstrated that foreground contribution is neglegible and 
therefore we decided not to add any foreground templates, 
CMB sky is gen erated according to the W MAP best-fit 
power spectrum jBennett et alJl2003l) and convolved with a 
symmetric Gaussian beam with the BEAST angular resolu- 
tion; 30' in Ka ba nd and 23' in Q band. Maps are pixelized in 
HEALPix format llGorski et al, 1999) with Nm^ = 512, Ob- 
serving these maps, we produced Time Ordered Data (TOD) 
for each BEAST channel from which we created maps in 
Ka and Q bands following the same reduct i on pro cessing of 
the real BEAST data (see iMeinhold et alJ (|200l) for map- 
making process details). 

As for instrumental noise simulation, we did not pro- 
duce noise time streams for each of BEAST detectors but 
we adopted a different recipe in order to have noise maps 
with the same statistical properties as the actual data. 

We generated white noise realizations with the same 
rms per pixel of BEAST maps in the two bands. It is clear 
that this is not enough due to the non-negligible level of 
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Table 2. Simulated maps S/N ratios 



S/N 


Ka 


Q 


30' 


~ 0.09 


~ 0.57 


40' 


~ 0.28 


~ 0.81 


60' 


~ 0.33 


~ 0.85 



Table 3. Number of "good" CMB reconstructions from 100 sim- 
ulated maps for each angular scale using p, g and t function in 
FastICA algorithm. See text for explanation. 



1// noise fe.o. iMeinhold et al.ll200^ which makes the noise 
clearly not -white. We made use o f the "difference" maps as 
derived bv IMeinhold et all (120051) to extract the noise angu- 
lar distribution. We then expanded in spherical harmonics 
both the white noise and the "difference" maps obtaining 
the harmonic coefficients and af^ respectively. These 
coefficients are combined according to: 



(3) 



and we then generated noise maps with these new coef- 
ficients. In this way noise simulated maps have also the same 
angular power spectra as the actual processed BEAST maps. 
We repeated this procedure for each of the two BEAST 
frequency bands. Finally, after smoothing to 30', 40' and 
60' angular resolution, we added CMB and noise simu- 
lated maps together at each frequency, obtaining simulated 
BEAST maps. Sub-pixel noise effects are negligible due to 
the high angular resolution, and the smoothing reduces any 
possible residual effect, possible In Table |21 we report the 
S/N ratios of the simulated maps. We underline that sim- 
ulated Q ma ps smoothed to 30' ha ve the same S /N ratio 
estimated bv IMeinhold et alJ for the BEAST 30' Q 

map derived with full processing of the data. This fact is 
a direct verification of the success of our recipe for noise 
simulations. 

Subsequently we applied FastICA to the simulated 
maps, after removing the region with |&| ^ 17.5° as done 
for the actual data. For every run we derived correlation co- 
efficients of both ICA maps with the input CMB. Therefore 
the ICA map of the two with the higher correlation coeffi- 
cient is a possible CMB reconstruction. 

Figure |2| shows the relation between correlation coef- 
ficients and frequency scaling of the possible CMB recon- 
struction at the three different angular resolutions for the 
three non-quadratic forms assumed by ICA. It is interest- 
ing to note that in corrispondence of the expected CMB 
frequency scaling (1.022), we observe the higher correlation 
coefficients and then the best recovered CMB. We therefore 
use frequency scaling as a figure of merit for the CMB recon- 
struction. Furthermore we observe from Fig.j^jthat increas- 
ing angular resolution, increases also values of the correla- 
tion: CMB reconstruction becomes better. This is expected 
because the S/N ratio increase with angular resolution (see 
Table HJ. 

Finally we used this relationship between correlation 
and frequency scaling to establish which CMB reconstruc- 
tions are reliable. For every resolution we selected the mini- 
mum correlation value that characterizes a good reconstruc- 
tion. This correlation coefficient corresponds to a maxi- 
mum value of frequency scaling s. In Tablej^we report the 
chosen values of correlation coefficient and frequency scaling 
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(selecting regions near the peaks in Fig. |5J and the number 
of "good" recove red CMB maps. As alre ady noted in pre- 
vious works (e.g. iMaino et al"ll200l l2003l) FastICA results 
with t function are the worst, while g performs better in the 
astrophysical context. Increasing angular scale increases the 
number of "good" CMB reconstructions and decreases the 
differences between p, g and t. 

Looking at the CMB recovered from the data, we verify 
that the frequency scalings reported in TableQlie within the 
range of values that identify at every resolution a "good" 
reconstruction, with the exception of t results at 40' and 
60'. This is a further indication of the fact that FastICA is 
able to extract reliable CMB signal. 



4 NORMALISATION OF THE CMB 
COMPONENT 

In general FastICA recovers a copy of the original signal 
i.e. it is not able, in principle, to recover the variance of 
the underlying sources. Therefore we have to normalize the 
CMB component recovered from BEAST data. Generally we 
can derive the right normalisation factor directly from the 
FastICA outputs, but in this case, due to the poor S/N 
ratio, we must again use our Monte-Carlo simulations. In 
fact for each CMB recovered from simulations the scale fac- 
tor is just the ratio between reconstructed and input CMB 
maps rms. The output CMB map is in Ka band antenna 
temperature and we compare it with the simulated CMB 
map in Ka band. However noise in the reconstructed CMB 
is quite important and we have to estimate and subtract it 
(at least in terms of rms) . Although the FastICA algorithm 
is higly non-linear, the data model is linear i.e. sources are 
obtained with a linear combination of the input data. In this 
way we can quite easily estimate the noise contribution in 
the reconstructed components by exploiting the separation 
matrix W elements pertinent to the CMB component. We 
then subtracted, for each simulation, noise rms from the the 
output CMB rms to obtain the exact recovered CMB rms 
to compare with the input CMB rms. 

Results are shown in Fig. |5] where normalisation factor 
is reported as function of the frequency scaling at the three 
different angular resolutions. There is a clear relationship 
between normalisation factor and frequency scaling. Indeed 
this relation is almost linear within a frequency scaling Smax 
and a normalisation factor Nmax ■ Such values decrease when 
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Figure 3. Normalisation factor as a function of tlie frequency scaling for 30', 40' and 60'. Line-dot lines indicate normalisation factor 
equal to 1 and the expected frequency scaling. The bigger symbols show the interpolation with the frequency scalings of the CMB 
recovered from BEAST data, cutting the galactic plane for |6| ^ 17.5° and using p (asterisks), g (diamonds) and t (triangles) function. 
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Figure 4. Normalisation factor as a function of frequency scaling 
for 30' (asterisks), 40' (diamonds) and 60' (triangles). 



5 POWER SPECTRUM 

We estimated the angular power spectrum of the FastICA 
CMB compone nt from BEAST data and compared it with 
that derived bv lO'Dwver et al] (|200^ for the analysis of Q 
band BEAST data. We extracted the spectrum choosing an 
angular resolution of 40'. This is a good compromise be- 
tween S/N ratio and signal level which is affected by the 
10 Hz high-pass filter applied to the data. We considered 
the more conservative galactic cut (|6| ^ 17.5°) and used 
the g function in the ICA algorithm which has the better 
recovered frequency scaling (see Table 0. 

To extract the CMB power spectra we adopted the 
MASTER method (H ivon et al. 2002) that was also used by 
lO'Dwver et alj J2005 l). MASTER returns a binned pseudo- 
Ce estimator allowing for de-biasing the power spectrum for 
the effects specific of the experimental CMB observation, 
such as sky-cut, scanning strategy, data processing and in- 
strumental noise. This is expressed by the following data 
model: 



angular scale increases since FastICA performs better (see 
TableOJ and scalings identifying a "good" CMB reconstruc- 
tion are within Smax- Furthermore when the CMB compo- 
nent has the expected frequency scaling, it has also the right 
normalisation, with scale factor equal to 1. 

In Figure |3 we also report the frequency scalings ob- 
tained from BEAST data at the three angular resolutions 
and for the three non-quadratic functions. We derived the 
normalisation factor by simple linear interpolation of this 
relation in the points corresponding to the actual frequency 
scalings. 

Finally we observe that this relation does not depend 
on the functions p, g and t, and not even on angular reso- 
lution for scaling smaller than the theoretical one as shown 
in Fig.Bl 



Ci = Y^ Mw Ft, Bl Cf + {Ne), (4) 

where Ci is the observed power spectrum, C}^ is the the- 
oretical one. The i3| term includes both instrumental and 
pixel window functions and the kernel Aiui accounts for the 
mode-mode coupling between different modes due to the in- 
complete sky coverage and depends on the actual shape of 
the observed sky region (so it can be computed once for 
all). The other terms are calibrated against Monte Carlo 
simulations. In particular with simulations of CMB only ob- 
servations we compute the instrumental transfer function 
Fi, which accounts for data processing effects. Instrumental 
noise only simulations are needed to estimate the average 
noise angular power spectrum (A'^^), while from simulated 
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skies (CMB + noise) we derive errors on o ur final power 
spectrum estimation (see lHivon et alJ i2002l) for details). 

For a proper application of MASTER to a FastICA 
CMB map, we must take into account that an ICA CMB 
map is indeed a linear combination of two data maps: 
the one in Ka-band and the one in Q-band, i.e. ycMB = 
WKa^Ka + WQ'x.Q, where the weights WKa and wq are de- 
rived from the ICA separation matrix W. In general, due 
to the different S/N ratios, we have obtained wq > WKa 
(for example in the chosen case we have WKa ~ 0.20 and 
Wq ~ 1.04), however the ICA CMB power spectrum is af- 
fected by experimental observation effects in both the bands. 
In order to evaluate such effects we made use of the Monte- 
Carlo simulations already performed at 40'. First we used 
simulated Ka-band and Q-band CMB to estimate the instru- 
mental transfer functions for Ka and Q bands separately, 
-F/^" and F'^ respectively. We then computed the final ICA 
transfer function as Fi = w^^fI^"" -f wqF^ . Intrumental 
noise, and its angular power spectrum, in the reconstructed 
CMB component is obtained in a similar manner: we use the 
same coefflcients WKa and wq to properly combine the Ka- 
band and Q-band instrumental noise realizations. Finally 
the same weights are used in signal plus noise simulations in 
order to derive final error on the MASTER power spectrum. 

We extracted the binned power spectrum choos ing the 
same multipole bins used bv lO'Dwver et al] ^2005), with 
A£ — 55. We estimated the ICA power spectrum for nmlti- 
poles ^ 400 since signal is significantly suppressed at higher 

multipoles due to the selected 40' smoothing. 

Finally, since FastICA does not recover a CMB signal 
with the right variance, we normalised the spectrum using 
the scale factor derived from interpolation of the "normal- 
isation factor - frequency scaling" relation, as described in 
the previous Section (see Fig. ^ . 

The resulting 40' FastICA CMB power spectrum is 
shown in Fig. |3 compa red with the 23 ' Q-ba nd map power 
spectrum estimated bv lO'Dwver et ahl l|200^ and with the 
best-fit WMAP model. The agreement between the spectra 
is good. In particular the two BEAST spectra agree within 
1 — 0". In Table 2] we report the d and associated 1 — a 
errors for the two BEAST power spectra. Furthermore, both 
from Fig.lHland Tablegl we can observe that ICA spectrum 
error bars are smaller than those of Q-band power spectrum. 
This is due to the smoothing at 40' , which reduces the noise 
contribution. 

The spectra agreement is a strong indication of the suit- 
ability of the FastICA CMB signal, and also of the goodness 
of the adopted normalisation procedure. 



6 CRITICAL DISCUSSION AND 
CONCLUSION 

In this paper we applied FastICA algorithm to real CMB 
data from the BEAST experiment. This is a ground-based 
experiment operating from the UC White Mountain Re- 
search Station (CA) at an altitude of 3.8 km that produced 
partial sky maps in two frequency bands (Ka and Q) with 
angular resolution of 30' and 23' respectively. 

One of the FastICA requirements is that the instru- 
mental noise has to be Gaussian and uniformely distributed 
on the sky. This is not the case for BEAST, which clearly 



shows 1/f noise and non-uniform integration time due to the 
observing strategy. The 1/f noise has been accounted for by 
applying a high-pass filter to the Time Ordered Data in the 
map-making process. This of course alleviate the impact of 
non-white noise but also reduces sky signal on large angular 
scales. Furthermore the si gnal-to-noise rat i o, as estimated 
from "difference" maps bv lMeinhold et al] i2005l) . is quite 
poor being ~ 0.11 for Q-band and even lower in Ka-band. 

Another limitation for FastICA applicability is that 
different frequency channels have to be at the same angular 
resolution. This forces us to further convolve our data set. 
We choose three different values for resolution: 30', 40' and 
60'. We did not apply a more aggressive smoothing since 
the high-pass filter effect on data is a clear suppression of 
signal on larger scales. Smoothing data allows us to reach a 
slightly better S/N ratio which helps in the application of 
FastICA. 

All these constraints have the consequence that Fas- 
tICA always extracts from BEAST data, one physical com- 
ponent while the other is clearly noise related. Furthermore 
in order to extract a CMB component we have to cut out 
the galactic plane where galactic emission is dominating over 
CMB since otherwise this would prevent us from properly 
reconstructing CMB. Nevertheless, after galactic cut, Fas- 
tICA recovers a CMB-like component, but with frequency 
scalings (Table quite different from the theoretical one. 
This is again due to the relatively high instrumental noise, 
which alters CMB reconstruction. Despite this first bad indi- 
cation, with further analysis we verified that this component 
is indeed dominated by CMB anisotropics. 

In order to test our CMB results quality we ran Fas- 
tICA on 100 Monte Carlo simulations of Ka- and Q-band 
data at the three selected angular resolutions. These have 
been simulated by creating fake CMB skies with angu- 
lar power spectrum from the best-fit model from WMAP 
jBennett et alJ l2003^ . observing these skies according to 
BEAST scanning strategy and reducing data with the same 
pipeline applied to real data. We finally superimposed in- 
strumental noise with the same statistical and spatial prop- 
erties as the actual data. We derived correlation coefficients 
between CMB ICA maps with the input CMB and stud- 
ied the relation between these coefficients and the recovered 
scaling frequency, in order to use scaling as a figure of merit. 
For every resolution we selected the maximum scaling al- 
lowed for a "good" CMB reconstruction and the comparison 
of these values with the results out of BEAST data confirms 
that FastICA indeed recovers a reliable CMB component. 
Furthermore the relation shows the increasing reconstruc- 
tion quality with angular scale, since the increasing S /N ra- 
tio, and the different FastICA performance with the three 
non-quadratic functions, with the better and the worse re- 
sults obtained using g and t functions respectively. This final 
indic ation agrees with previous works (Maino ct al. 20o3) 
E)Ol . 

Since FastICA is not able to recover the variance of 
the independent components (it recoves a "copy" of the in- 
dependent underlying components), we again used Monte 
Carlo simulations to obtain a normalisation procedure for 
the CMB component. In fact in this case the scale factor is 
just the ratio between output and input CMB rms. For each 
resolution we found a clear relation between scale factor and 
frequency scaling that is almost linear within certain values 
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Figure 5. CMB power spe ctra from BEAST dat a: stars show the spectrum extracted from FastICA 40' CMB map, diamonds show the 
23' Q band map spectrum iO'Dwver et aLl2^^,'^^ . ICA spectrum is shifted by A£ = 4 for clarity. Solid line is the best-fit W MAP power 
spectrum. 



Table 4. The BEAST Cg and la error values (in AtK^). 



bin 


FastICA 


map 


Q map 






i(£ + l)Ce/2n 


Icr error 


£(£+ l)Q/27r 


1(T error 


139-193 


3865 


±425 


3776 


±552 


194-248 


4628 


±497 


4744 


±781 


249-303 


3783 


±535 


3597 


±782 


304-358 


4237 


±528 


3374 


±625 


359-413 


2997 


±703 


1829 


±969 



of scaling and normalisation factor. The decreasing of such 
values with angular scale indicates again the corresponding 
improvement of FastICA performance due to the better 
S/N ratio. Furthermore the relation does not depend on the 
non-quadratic function and shows that those CMB recon- 
structions with the expected frequency scaling have also the 
correct normalisation {e.g. equal to 1). The normalisation 
factors for FastICA results out of BEAST data are derived 
by interpolation of this relation at the derived frequency 
scalings. 

Finally we extracted the FastICA CMB a ngular power 
spectrum adopting a MASTER approach iHivon et alJ 
[2002) and normalised it with the proper scale factor. We 
found a very good agreement with our results and the best- 
fit WMAP model and also wit h the spectrum estim ated 
from the BEAST Q-band map dO'Dwver et al.ll2005l) . al- 
though on a limited multipole range because of the extra- 
smoothing applied to the data. This spectra agreement con- 
firms the rehability of the CMB extracted by FastICA and 
validates our normalisation procedure. 

Our analysis, toge t her w ith that on DMR data per- 
formed bv lMaino et al.l ll2003ll . confirms the very good per- 
formance of blind algorithms like FastICA in extracting a 
CMB component even from noisy data on a small patch of 
the sky like BEAST ones. Therefore we think that blind 



algorithms are valid tools for present and future CMB ex- 
periments providing information on the independent compo- 
nent in the actual observed sky signal which could be used to 
feed much more complex algorithm like Maximum Entropy 
Method. This is particularly relevant for future CMB polar- 
isation experiments where we will be forced to work with 
low S/N ratios and wher e our knowledge of polarisation for 
foregrounds is still poor llStivoh et alJl2005^ . 
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